On the transfer of entanglement from a two-mode squeezed state to a pair of qubits 
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There have recently been interests in transferring entanglement between two quantum systems in 
different Hilbert spaces. In particular, the study of entanglement transfer from a continuous- variable 
to a qubit system has a primary importance due to practical implications. A continuous-variable 
system easily propagates entanglement while a qubit system is easy to manipulate. We investigate 
conditions to entangle two two-level atoms using a broad-band two-mode squeezed field driving the 
cavities where the atoms are. 
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PACS numbers: 42.50.p, 03.65.Ud, 03.67.-a, 42.50. Pq 

Introduction - The possibility to exploit the correla- 
tions transferred from a two-mode squeezed state of light 
to two distant atoms, initially prepared in a separable 
state, appeals for the purposes of quantum communica- 
tion and distributed quantum computation Q]. On the 
other hand, some of the aspects that characterize the dy- 
namics of two-level atoms in a squeezed reservoir have 
been studied and exploited for quantum state engineer- 
ing 0. Very recently, Kraus and Cirac have sug- 
gested a cavity-quantum-elctrodynamics (CQED) system 
in which two trapped atoms in respective optical cavities 
interact with an external broad-band two-mode squeezed 
state. The scheme has been exploited to produce an en- 
tangled atomic state. 

Starting from this system, we address the condi- 
tions that have to be satisfied to perform effective 
entanglement-transfer from a continuous- variable system 
to a pair of qubits. This process is intriguing because 
it combines the manipulability of a qubit system to the 
possibility of entanglement propagation via a correlated 
state of light. We analyze the entangling properties of 
the master equation that describes the dynamics of the 
system in [jjj. We find bounds on the degree of purity 
and entanglement of the two-mode driving field, assess- 
ing the corresponding effects on the reduced state of the 
two atoms, both for dynamical and steady conditions. 
We consider the cavity decay rate and the atomic spon- 
taneous emission for the efficiency of the scheme. 

Master equation - Here, we briefly review the model 
recently proposed in 3]. Two identical two-level atoms 
are respectively trapped in spatially separated cavi- 
ties, which are externally driven by a broad-band two- 
mode squeezed state of light characterized by the phase- 
insensitive factor TV and the phase-sensitive one M. The 
two identical cavities have an energy decay rate k. Un- 
der the assumption that the bandwidth Auj ext of the 
squeezed field is greater than the cavity decay rate k, 
the Born-Markov approximation is applied. The master 
equation reads (H = 1) 



d t p = -i[H al + H b2 , p] + t c 



(1) 



with H a i = fi (pi o, + a) a 1 ) a Jaynes-Cummings inter- 
action Hamiltonian between atom 1 and cavity mode a 



Q, is the Rabi frequency. The Liouville operator C cav de- 
scribes the interaction of the cavity-mode fields with the 
external broad-band squeezed field as 



£cavP = K {(N + l)(apa^ — a^ap) + N{6^ pa — do^p)} 

ct=a 

+ 2nM(apb + bpa — bap — pivw) + h.c., 

where h. c. denotes hermitian conjugate of all the terms 
and p contains the degrees of freedom of both the atoms 
and the cavity fields. The contributions from the sponta- 
neous emission of the atoms from the excited state \e) ab 
have been neglected. However, they will be included 
later. 

For the case of a minimum-uncertainty squeezed state 
with a real squeezing parameter r, it is N = sinh 2 r, M = 

Using the squeez- 

abj | (assuming r 

real) Q it is possible to transform p as p = S(r)pS^ (r), 
then Eq. (JTJ becomes dtp = (£ a + £cav,o) P- Here we 

have defined L a p = — i[H a i + Ti-b2,p], where H a i = 
S(r)H a iS'(r) (analogously for Hb2) and, moreover, 
Ccavfl '■= C cav when N = M = 0. 

Trap state - It is straightforward to see that the Hamil- 
tonian H a i+Ti.b2 bas the eigenstate |*) = \\jj) 12 ® |00) a6 , 



sinhrcoshr and M 2 = N(N + 1) 
ing operator S(r) = exp 



N 



From 



(analogously for H^)- Here, <rj~ = (& 1 



) 1 (g\ and 



where \ip} 12 = (y 7 \gg) + y^jv+i 
the Liouvillian C cav ,o it is possible to identify the jump 
operators Ci = \/kci, C2 = \ftib that are responsible 
for the open dynamics of the cavity modes. Obviously, 

Tr (ELi & Ci) = as Q |*> = 0, so that |*) is 

a steady state solution with the jump rate equal to zero. 
Hence, we call |W) a trap state of the system. For a closed 
dynamics (k = 0), the only way to populate the trap state 
is to initially prepare the system in the trap state itself. 
In the limit of k < fl, on the other hand, the probability 
that the system, after its transient period, will end up 
in the trap state |^) is not negligible only if the system 
is actually prepared in this state. In general, we find 
that the population of the trap state changes according 
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to $«*|p|tf)) = E i= i, 2 (*|6<P<5?|*) oc k. The rate at 
which the trapping state is populated because of the open 
dynamics of the system is the decay rate k. Indeed, k is 
the rate at which the quantum correlations between the 
cavities can be built up. It is, thus, convenient to choose 
a value k ^> f2 to increase the probability that the atomic 
steady state is populated. This condition corresponds to 
a bad-cavity regime. However, k can not be increased at 
will because the relation Aio ext ^ « represents, in gen- 
eral, a strong constraint for the validity of Eq. 1(1)1. Typi- 
cally, in fact, it is Au> ex t — k/6 ~ 2nx 12 MHz, in experi- 
ments in which squeezed light is coupled to a cavity-atom 
system [f|. However, an analysis in the finite-bandwidth 
case is out of the scope of the current work and the con- 
dition Au) ex t 3> K will be assumed throughout this work. 

Bad- cavity limit - We now assume the bad-cavity limit 
k 3> f2, I\ where T is the atomic spontaneous emis- 
sion rate. Values (fi, T) /2n ~ (20, 3.5) MHz allow for 
the squeezed state to build up inside the cavity p|. In 
this weak-coupling regime, the dynamics of the cavity 
modes interacting with the external driving fields is much 
faster than their interaction with the atoms. In these 
conditions, the atoms see the cavity modes in a steady 
state that is not affected by the atom-cavity dynamics. 
Thus, for the interaction time such that r = nt 3> 1, 
we take the steady state of the cavity modes to be a 
two-mode squeezed state as in the case without atoms 
being in the cavities. We are interested in the atomic 
evolution so to eliminate the cavity modes. In virtue 
of the above consideration, the adiabatic elimination of 
the cavity fields is possible by defining the projection 
operator as Vp : p — > Vp = p ss ® Tr ab (p) = 
Pss ® Pi2- Here pi2 is the atomic density operator. It 
is easy to prove that VC a Vp = and it leaves the 
master equation in the adiabatically eliminated form 

d t pi2 = Tr ab {£ a / °° •»'£„ (p ss <g> p 12 ) dt}. To fur- 

ther reduce this equation we have to invoke the first 
Born-Markov approximation and exploit the relation 
£cav,o{Pss) = 0. Going back to the initial unsqueezed pic- 
ture, the adiabatically eliminated master equation takes 
the form 

2 

dtpn = - [{N + 1) (d--p 12 &+ - cr+a7p 12 ) 

i=l 

+<7~p 12 <7f - frz&iPu - pi 2 d-+a+) + h. c. 

(3) 

We have introduced the effective energy decay rate 7 = 
2Q 2 /k. This master equation is in agreement with the 
result found in Q . Eq. describes the effective decay of 
a pair of two- level atoms in a squeezed vacuum 0] ■ The 
same equation can be obtained relaxing the condition of 
a minimum-uncertainty two-mode squeezed driving field, 
even if the calculations involved are more complicated. 

In this analysis we can include the Liouvillian terms 
related to the scattering of the atoms. These terms 



depend exclusively on the atomic operators according 
to C sp ont,iP = r/2 (cr^pa+ - frr af p - pcr+a^ + h.c), 
with i = 1,2. The application of the projector opera- 
tor V to these terms of the master equation simply gives 
p ss <g) C sp ont,i(p 12)- Following the lines sketched above, 
we get a generalized master equation of the same form 
as Eq. © but with the parameters 7' = £1 2 (2 + T)/k, 
N' = CN/(1 + C) and M 1 = CM /{I + C) instead of 
7, N and M , respectively. Here we have introduced the 
cooperativity parameter C — 2Q 2 /Tk jij. With these new 
parameters it is always M' < \J N'(N' + 1), even for a 
minimum-uncertainty two-mode driving field 0- Thus, 
inside the cavity a pure two-mode squeezed state can 
never be obtained, in a realistic situation. 

Bound conditions for entanglement-transfer - It is 
widely believed that an interaction mediated by a non- 
Markovian reservoir can set entanglement between two 
subsystems because of the memory effects of the bath. 
However, recently it was shown that even a Markovian 
reduced dynamics is able to set quantum correlation be- 
tween two initially non interacting systems • Under the 
hypothesis of weak-coupling and Markovian nature, the 
reduced dynamics of two systems, A and B, interacting 
with a common bath can be cast into a master equation 
form as d t pAB = -i[H total, Pab] + CpAB, with 

CpAB = Da P (OaPABOp - - \OpOu, PAS I J • 

(4) 

Here, O a = a a <g> 1 for a — 1, 2, 3 and O a = l® er a _ 3 for 
a = 4, 5, 6, with 01,2,3 the Pauli matrices and D a 6 x 6 
matrix decomposed as 

D = (£ *) , (5) 

where A = A^,C = and B are 3x3 matrices. 
Following ref. 7], it is possible to characterize the en- 
tanglement capabilities of the bath-mediated interaction 
between the two atoms of the cavity-QED system we 
are treating here. Namely, the condition to entangle 
the atomic subsystem using the Markovian dynamics de- 
scribed by Eq. © is (u*Au) (v*C T v) < |u*fte(B)v| 2 , 
where u = (cos 29, — i, 0) T and v = (cos 2ip, i, 0) T carry 
information about the generic initial states of atoms 1 
and 2, unitarily rotated by the angles 6 and ip around 
the z axis of the Bloch sphere. For the case of Eq. 0, 
we identify 

/m_±i 1 \ /if 0\ 

A = C = 7 -| 2iYl+i B = 7 ■ 

V 00/ V 0/ 

(6) 

While these matrices depend solely on the form of the 
reduced master equation for the atomic subsystem, the 
condition to entangle the two-level atoms depends on the 
initially prepared atomic state via the vectors u and v, 
that carry information about 9, ip. If the initial state 
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of the atoms is taken to be \gg) 12 (gg\, the condition be- 
comes 

N 2 < M 2 (7) 

to entangle two atoms in cavities. Eq. Q has a clear 
physical interpretation: it is the sufficient and necessary 
condition for the entanglement of the driving field [10( . If 
we start with the atoms in their ground states, the entan- 
glement of the driving field represents, thus, the irrcmis- 
sible condition for the creation of entanglement between 
the atoms, despite the atomic spontaneous emission |l l| . 
To confirm the validity of this relation, we have solved 
the master equation ©, looking for the set of Bloch 
equations for the matrix elements pijhk — (ij\pi2\hk) 
(i,j,h,k = e,g) under the assumption M 2 <N(N + 1) 
and r = 0, even though including V ^ is straightfor- 
ward. Using the basis {|ee) , \eg) , \ge) , \gg)}i2, we find 
the following set of coupled, linear differential equations 

Peeee = 7 [ 2n^p eeee + Tl^(p e g e g + Pgege) ~\~ 2M Peegg\ j 



Pegeg = 7 ^ll 1 
Pgege = 7 



Pgege) \ Peeee 
Pegeg) H - Peeee 



Pe 



where, 



-7 [nlp ee gg - M(l - 2p, 



n 3Pegeg 
n 3P9ege 
ege ^Pegeg 



2M Peegg ] 
2Mp eegg ] 

)] 



(kN + I) and, by hermiticity, p e 



(8) 



Pgge 



All the other matrix elements are decoupled from these 
equations. The condition Tr^ipu) = 1 gives the equa- 
tion fOr Pgggg. 

Solving Eqs. iJHJ, we can find the dynamics of entan- 
glement between atoms 1 and 2. We use the entangle- 
ment measure based on negativity of partial transposi- 
tion (NPT), defined by E NPT {t, N, M) = -2\r(t, N, M), 
where A~ (t, N, M) is the negative eigenvalue of the par- 
tially transposed density matrix pji (T 2 indicates partial 
transposition with respect to 2) [8j]. NPT is a necessary 
and sufficient condition for entanglement of any bipartite 
qubit system . The results are shown in Fig. ^ 
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FIG. 1: Entanglement of pn after the reduced Markovian 
dynamics described by the master equation @. The entan- 
glement is plotted as a function of the dimensionless interac- 
tion time r — ■yt for N = 0.7 and three different values of M: 
M = 0.79 (solid curve), M = B as (0.7) = 0.902 (dashed curve) 
and M = ^N(N + 1) = 1.09 (dot-dashed curve). In (a), the 
initial state is \gg} 12 {gg\ while in (b), it is |ee) 12 (ee|. 

By inspection of Fig. ^a), it is apparent that, even 



for N < M, the long time behavior of the entanglement 
function can lead to a separable atomic state. This is 
due to the fact that the condition found does not give 
information about the steady-state entanglement. The 
criterion for the entanglement due to the interaction with 
a Markovian bath is, indeed, relative to the creation of 
quantum correlation in an initially separable state, with 
a positive gradient of £jfpr(t, N, M). In order to be sure 
that the Markovian reduced dynamics will lead to the 
creation of entanglement, the trend of the entanglement 
function at the initial interaction time t = has to be 
positive. This is the case for pi2(0) = \gg)i 2 {gg\ but not 
for pi2(0) = |ee) 12 (ee|, for example. For this initial state, 
it is d t £NPT(0, N, M) < 0, as can be seen in Fig. ^b), 
and the entangling condition leads to N + 1 < M, t hat is 
physically meaningless since N + 1 > y/N{N + l), V N. 
An analogous result is found for the asymmetric case in 
which one atom is prepared in \g) while the other is in 
|e). However, as stated before, this does not exclude the 
possibility that the atomic state becomes entangled later 
on. 

It is, thus, interesting to investigate the conditions that 
have to be fulfilled for the atomic steady state to be en- 
tangled. In particular, we are interested in finding the 
boundary value B SS (N) of the phase-sensitive parame- 
ter M beyond which we are sure that the atomic steady 
state is inseparable. To find B SS (N) we have to look 
at the asymptotic behavior of the entanglement function 
that turns out to be a non-decreasing function of the pa- 
rameter M, for a fixed value of N, as can be seen solving 
Eqs. JHJ and looking for the steady solutions. Thus, the 
condition 

Hm £ N pr{t, N, M)\ M=Bss{N) = (9) 

fully characterizes the boundary value for M that we are 
looking for. Using this, we find 



B SS (N) = -a+ y/a + N{N + l) 



(10) 



with a = l/4(Air) 2 and (Ax) 2 = N +1/2 the variance 
of the in-phase quadrature of the bath. It is B SS (N) < 
y/N(N + 1), V7V. The boundary, Eq. {THJ, is indepen- 
dent from the initial preparation of the atoms. Thus, for 
the CQED system considered in this work, the atomic 
steady state is entangled when the two-mode driving field 
satisfies the condition B ss (N) < M < ^/W{N+~T). Out- 
side this range, the steady state is separable and the 
atomic subsystem may be entangled depending on the 
initial preparation of the atoms. This is one of the main 
results of this paper. 

Eq. I|10|) represents a very strict condition on the prop- 
erties of the driving field. The experimentally avail- 
able source of squeezed light are, indeed, quite bright 
(that means large values of N) and it is easy to see that 
B SS (N > 1) ~ y/N(N + 1), dramatically shrinking the 
range of values of M in which the atomic steady state is 
entangled. 
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FIG. 2: In (a), S L (t,N,M) is plotted against r = jt for 
N = 0.7 and M = y / N(N + 1) = 1.09 (solid line), M = 
1.05 (dashed line) and M = 0.8 (dot-dashed line). In (b), 
S L (oo,N,M) as a function of M for N = 0.7 (dashed line) 
and N = 0.9 (solid line). 

Finally, it is worth stressing that the the range of val- 
ues of M in which the steady state is pure reduces to 
the point M = y/N(N + 1). We look for the boundary 
value at which the atomic steady state becomes pure. 
As a measure for purity we take the linearized entropy, 
that for a two-qubit system is defined by Sl(N, M, t) = 
4/3 (l — Tr 12 {P12}) and ranges from (pure states) to 
1 (maximally mixed ones). Only the interaction with 
a pure squeezed reservoir realizes a pure atomic steady 
state. These results, for /012(C)) = |<?g) 12 (<?sl, are shown 
in Fig. [21 In (a) Sl is plotted versus the dimensionless 
interaction time for the case of a pure squeezed state 
(solid line), for M = y/N(N + 0.9) (dashed line) and 
M = y/N(N + 0.2) (dot-dashed line). In (b), we show 
the behavior of the linear entropy as a function of the 
phase-sensitive parameter M for N — 0.7 (dashed line) 
and N = 0.9 (solid line). It is apparent that the higher 
is N, the more mixed is the atomic state, for a fixed M. 
This is because higher values of N correspond to an in- 
crease of the thermal character of each field mode, in a 



two-mode squeezed state. On the other hand, increasing 
the squeezing properties of the driving (i.e., increasing 
M) leads to a rapid decrease of the mixedness of the state. 
This is because the populations p e geg = Pgege decay if we 
increase M, a signature of the two-photon correlations in 
a highly squeezed two-mode state Q- 

Remarks - In this paper, we investigated about the 
mechanisms for the transfer of entanglement from a 
continuous-variable state to a qubit system. Starting 
from the recent proposal Q , we have shown that the ef- 
ficiency of the scheme depends on the cavity decay rate. 
We have addressed the conditions for the entanglement of 
an atomic subsystem induced by the Markovian interac- 
tion with a broad-band two-mode squeezed bath, both in 
the dynamical and steady conditions. We have found the 
sufficient and necessary conditions for the driving field to 
entangle two atoms in cavities. The entanglement of the 
driving field is only a necessary condition to entangle the 
two atoms, in steady state. We found the boundary value 
of M beyond which the atomic steady state is entangled 
and pure. The effect of the atomic spontaneous emission 
rate T can be included in our analysis and Eq. (|10() is 
still valid replacing N with TV'. Our analysis gives some 
insight in the role played by purity and quantum corre- 
lation of the squeezed state in an entanglement-transfer 
process. 
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